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ABSTRACT 

We present Gemini Multiobject Spectrograph integral field spectroscopy of the extended emission-hne region 
associated with quasar 3C 249.1. The kinematics of the ionized gas measured from the [Olll] A5007 line is 
rather complex and cannot be explained globally by a simple dynamical model, but some clouds can be modeled 
individually as having locally linear velocity gradients. The temperatures of the ionized gas appear uniform 
(varying from ^^12000 to 15000 K), while the densities vary from a few tens to a few hundreds cm^^^. The 
emission mechanism of all of the emission clouds, as indicated by the line-ratio diagnostics, is consistent both 
with "shock + precursor" and pure photoionization models. The total mass of the ionized gas is on the order of 
10^ M0. We estimate the bulk kinetic energy and momentum of the extended emission-line region of 2.5 x 10^^ 
ergs and 10^° dyne s, and a dynamical timescale of ~10Myr. By comparing the injection rates of kinetic energy 
and momentum of different galactic wind models with the observation, we argue that the emission-line clouds 
are most likely a direct result from the feedback of the quasar We also discuss the nature of the extended X-ray 
emission surrounding the quasar. 

Subject headings: galaxies:interactions — quasars; individual(3C 249.1) 



1. INTRODUCTION 

About a third of low-redshift (z < 0.45) steep-spectrum 
radio-loud quasars are surrounded by luminous emission-line 
region s that often extend to several tens of kpc from the nu- 
cleus dStockton & MacKentvlfi^P ). The morphologies of 
these extended emission-hne regions (EELRs) are often com- 
plex and clearly display a non-equilibrium situation show- 
ing strings of knots and filaments straddling tens of kpc, and 
they are in general uncorrected with either the distribution of 
stars in the host galaxy or the structure of the extended radio 
source, if present. 

The origin of the gas comprising these EELRs remains 
uncertain. The apparent correlation between the occurrence 
of EELRs and evidence for strong interaction inspired the 
idea th at the gas is tidally disrupted material from a disk 
dStockt on & MacKentv 1987). But this is unlikely because 
( 1 ) there is generally no correlation between stellar tidal fea- 
tures and the distribution of the gas, and (2) some confine- 
ment mechanism is required for the gas to retain the den- 
sities implied by the emission-line spectra for more than 
a very brief time 10'*-10^ years; Fabian etal. 1987; 
IStockton et al. 2002). It has also been suggested that the 
gas comes from a cooling flow in a hot halo surrounding 
the quasar (Fabian et al. 1987). This scenario has been ruled 
out by deep Chandra X-ray observations of four quasars 
showing strong EELRs, since the hot (10^ K) halo gas from 
which the warm emission-line gas is suggested to condense 
is not detected, indicating strong cooling is not taking place 
([Stockton et al. 2006). And, in any case, the standard cooling- 
flow paradigm seems to have been ruled out by recent high- 
resolution X-ray spectroscopy (Peterson et al. 2003). 

If the gas cannot come from the intergalactic medium 

' Based on observations obtained at the Gemini Observatory, which is 
operated by the Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the Gemini partner- 
ship: the National Science Foundation (United States), the Particle Physics 
and Astronomy Research Council (United Kingdom), the National Research 
Council (Canada), CONICYT (Chile), the AustraUan Research Council (Aus- 
tralia), CNPq (Brazil) and CONICET (Argentina). 



(IGM) through a cooling flow, then it must originate from the 
interstellar medium (ISM) of the quasar host galaxy. Despite 
the non-detection of hot gaseous halos by the Chandra ob- 
servations, highly structured X-ray emission regions are seen 
around two of the four quasars surveyed. One of the plau- 
sible explanations for the X-ray emission is that it is ther- 
mal bremsstrahlung from high-speed shocks. On the other 
hand, detailed photoionization modeling indicates that a low- 
density phase ('^2 cm^"^) and a high-density phase (~500 
cm^-^) are needed to explain the optical spectrum of one of 
the EELR clouds ai'ound 4C 37.43 (El;|§tockton et al. 2002). 
The two phases cannot be in pressure equilibrium because 
their temperatures are almost the same. A reasonable pos- 
sibility is that the high densities in the gas are continuously 
regenerated by shocks propagating through the surrounding 
medium. 

The high-speed shocks could be driven by a galactic 
wind, which itse lf results from the feedback from the quasar 
(Di Matteo et al. 2005) or a vigorous starburst in the host 
galaxy (see Veilleux et al. 2005 for a review). This evidence 
of ongoing galactic winds combined with the morphological 
properties of EELRs thus suggest a scenario where the EELRs 
are a direct result of a superwind — the gas is originally the 
ISM in the host galaxy and has been blown out and shocked 
by a galactic wind; the gas is then photoionized by the UV 
continuum from the quasar and/or ionized by the shocks. 

3C 249.1 (z = 0.31) is a powerful quasar, also known as 
PG 1 10 0-1-772. Around the q uasar, distinct X-ray emission 
regions (Stockto n et alJl20 0^) and one of the most luminous 
optical EELRs at 2: < 0.5 (Stockton & MacKentv 1987) are 
seen extending to radial distances of tens of kpc, making it 
an ideal candidate for an in-depth study. In this paper, we 
take advantage of the simultaneous spatial and spectral cover- 
age of integral field spectroscopy (IFS) to further explore the 
origin of the EELR and the extended X-ray emission around 
3C 249.1. Throughout we assume a flat cosmological model 
with Ho ^ 70 km s^^, f7„ = 0.3, and f^A = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 
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3C 249. 1 was observed with the Integral Field Unit (IFU) of 
the Gemini Multiobject Spectrograph (GMOS) on the Gemini 
North telescope. The queue observations were executed in the 
early half nights of April 7 and 8, 2005 (UT). Three exposures 
of 2825 s were obtained with B600/G5303 grating at a central 
wavelength of 6537 A. Between exposures the telescope was 
offset by about 0'.'4 to improve the spatial sampling. The half- 
field (one slit) mode was used to cover the emission lines from 
[Ne V] A3426 to [O III] A5007. With this setting, we had a 
field-of-view (FOV) of 3'.'5 x 5", a wavelength range of 4400 
to 7200 A, a dispersion of 0.46 A pixel" ^ and an instrumental 
full width at half-maximum (FWHM) of 3.5 pixels (1.6 A). 
The seeing varied from 0'.'8 to I'.'O. 

The data were reduced using the Gemini IRAF data re- 
duction package (Version 1.8). The data reduction pipeline 
(GFREDUCE) consists of the following standard steps: bias 
subtraction, cosmic rays rejection, spectral extraction, flat- 
fielding, wavelength calibration, sky subtraction and flux cal- 
ibration. Spectra from different exposures were assembled 
and interpolated to construct individual datacubes (x,y. A) with 
a pixel size of 0'.'05 (GFCUBE). Differential atmosphere re- 
fraction was then corrected by shifting the image slices at each 
wavelength to keep the location of the quasar constant. The 
datacubes were then binned to 0'.'2 pixels, which is the origi- 
nal spatial sampling of the IFU fiber-lenslet system. 

Since our study focuses on the emission line gas, it is de- 
sirable to remove the light of the quasar from the datacubes. 
For 3C 249.1 the quasar component can be cleanly removed 
from the datacubes by assuming that all of the continuum in 
the spectra is scattered light from the quasar, because (1) pre- 
vious studies have shown that the host galaxy of 3C 249.1 is 
very faint in optical compared with the quasar itself, and (2) 
extended emission line clouds usually show essentially pure 
emission-line spectra {e.g., Stockton et al. 2002). The detailed 
procedure was as follows: The model quasar spectrum was 
derived by cubic-spline smoothing the average spectrum ex- 
tracted from an aperture centered on the quasar Then the ratio 
of the spectrum from each pixel to this template quasar spec- 
trum was fit with a smooth curve comprising a small num- 
ber of cubic spline pieces, and the continuum was removed 
by subtracting the quasar template normalized by this curve. 
For each spectrum, pixels near strong emission lines were ex- 
cluded from the sample of the continuum fitting, and the loca- 
tions of these spectral regions were determined by measuring 
the wavelength of the [O III] A5007 line from each residual 
spectrum. These steps were repeated a few times to achieve 
best residual datacubes. Finally, the residual datacubes were 
smoothed to a common spatial resolution, and merged to a 
final datacube. We generated exposure maps for each wave- 
length to account for the exposure time variation caused by 
the dithering and the differential atmosphere refraction. 

3. RESULTS 

3.1. Kinematics 

Kinematics of the ionized gas can be measured from strong 
emission lines. Since a single 47 min exposure is enough to 
acquire good S/N in the [O III] A5007 line region, we derived 
the velocity fields from the single datacube that showed the 
best spatial resolution. Figure[3compares the broad-band im- 
ages created from the datacube and the HST WFPC2 image 
along with the radio jet. The velocity field is shown in Fig- 
ure|3in three separate velocity bands from —460 to H-590 km 
s"^ (relative to the quasar narrow-Une emission) to separate 



different clouds that are present along the same line of sight. 
The FWHM measurements were corrected for the ^70 km 
s"^ instrumental resolution. Thanks to the high spectral res- 
olution of the data set, at least eight EELR clouds at different 
velocities are identified from the datacube. In the HST im- 
age, two of the EELR clouds are blended together since they 
project onto each other, and many are seen as faint knots be- 
cause of the shallow exposure (300 s) and the truncation at 
about H-375 km s^^ of the F656N filter 

The EELR exhibits rather complex global kinematics which 
cannot be explained by a simple dynamical model. The three 
most luminous EELR clouds in the velocity range of —230 
to H-150 km s"^ show a velocity gradient from N to S that 
aligns with the clouds. This gradient may be an artifact of a 
common disc rotation inherited from the host galaxy; but the 
high velocity clouds at u > 200 km s"^ show a velocity gra- 
dient from E to W, which is almost perpendicular to the pre- 
vious one. Also, the southern cloud in the —230 to H-150 km 
s"^ panels shows a velocity increase of ~100 km s"^ along 
NE-SW direction, which is consistent with the elongation ob- 
served in the HST image, thus it can be ascribed to a simple 
rotating disk (at another different orientation); however, the 
line widths also decrease along the same direction, suggest- 
ing two discrete clouds are present as hinted in the HST im- 
age. So overall the kinematics are disordered, but on the other 
hand the velocity field seems locally ordered. 

There are apparently two different species of emission-line 
clouds — one showing relatively broad lines (FWHM ^300 
km s~^), the other showing narrower lines (FWHM ^^100 km 
s"^). The relatively broad line regions are more luminous, 
have smaller radial velocities {\v\ < 200 km s"^), and are 
identified as spatially resolved blobs in the HST image. In 
contrast, the narrower line regions tend to have larger radial 
velocities {\v\ > 300 km s"^) and are faint and seen as un- 
resolved knots in the HST image. These knots feature an H/3 
luminosity of > 10''" erg s"^. The most straightforward in- 
terpretation is that the narrower line regions are less massive 
clouds, therefore less luminous and more easily accelerated 
to a higher velocity; they are reminiscent of the "bullets" fre- 
quently seen in Herbig-Haro objects and local galaxies ex- 
hibiting superwinds {e.g., Cecil et al. 2001) but are orders of 
magnitudes more energetic. 

We notice that in the —230 to H-150 km s"' panels there is 
a region of very low surface-brightness emission (east of the 
bright emission cloud close to the center), but showing very 
broad lines (FWHM > 600 km s"'). It cannot simply be a re- 
sult of blended emission lines scattered from the bright clouds 
nearby, since from the regions between the clouds (refer to the 
middle right panel of Fig.|2ji we can tell that blending can in- 
crease the line width up to only ^400 km s"'. The fact that 
the region is along the radio jet (see Fig.^J hints that the gas 
is disturbed by the radio jet. 

3.2. Electron density and temperature 

The luminosity-weighted average electron density can be 
determined from the ratio of the [Oil] AA3726, 3729 dou- 
blet, and the electron temperature can be determined from 
the [O III] A4363/(A4959+A5007) intensity ratio {e.g., Oster- 
brock 1989). The S/N ratio from a single 0'.'2 pixel is not 
high enough to perform this kind of measurement, so we com- 
bined spectra from individual emission regions using an aper- 
ture of 0'.'4 in radius and a centered PSF as a weighting func- 
tion to achieve optimal extraction. As the last step, we cor- 
rect the spectra for the line-of-sight Galactic extinction {Av 
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Fig. 1.— Le//— archived HS T/WFPC2 F656N i mage of 3C 249.1 (the redshifted [OlII] A5007 falls into this filter), overlaid with contours of a VLA image 
at 5 GHz and Qf.'SS resolution IGilbert et alj|200^ . Contours at 0.5 X (1, 2, 4, 8) mJy beam^^. The peak of the radio continuum has been registered to the 
position of the quasar. The two X-ray emission regions (xa & xb) are marked in red and labeled. The gray-scale images in this and the middle panel have 
been allowed to wrap around to show both low-surface-brightness detail and high-surface-brightness peaks. Middle & Right — GMOS/IFU [O 111] A5()07 radial 
velocity broad-band channel maps with contours of the HST image and contours of the radio jet overlaid, respectively. The radial velocity range (relative to that 
of the nuclear narrow line region, z = 0.3117; negative velocities are blue shifted) is shown in each panel. The rectangles deHneate the 3'.'5 X 5" FOV of 
GMOS/IFU. The crosses indicate the position of the quasar, which has been removed from the datacube (for details of the PSF removal, see § ??). 
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Fig. 2.— Velocity field of 3C 249.1 EELR derived from the [Olll] A5007 region in the GMOS/IFU datacube. The fiiree columns, from left to right, are line 
intensity, radial velocity (relative to that of the nuclear naiTow line region) and velocity width maps. To separate different clouds that are present along the same 
line of sight, the velocity field is shown in three separate velocity bands (—460 to —300 km s , —230 to -1-150 km s~^ and +180 to -1-590 km s~ from the 
top row to the bottom row). The rectangles delineate the 3'.'5 X 5 " FOV of GMOS/IFU. Pixels are 0'.'2 squares. The crosses indicate the position of the quasar, 
which has been removed from the datacube. For details see S I3.1I 



= 0.112; Schlegel et al. 1998), and for the intrinsic redden- 



ing by dust associated with the clouds, which is estimated by 
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matching the measured H7/H/? ratio with the value predicted 
by case-B recombination. I n both cases we used the standard 
Galactic extinction curve o f iCardelli etaflfTosl . Unlike the 
kinematic analysis, these line ratios were extracted from the 
merged datacube in order to obtain better S/N for faint lines. 
Figure |3] shows the extraction regions, and Fig. |3 shows the 
extracted spectrum of EELR-6 as an example. The spectra 
usually show either asymmetric line profiles or two compo- 
nents with different radial velocities. For both cases the line 
profiles usually can be fit quite well with two Gaussians. Fig- 
ure|5]shows an example of our two-Gaussian model fit to the 
line profiles. The spectrum shows two components at differ- 
ent radial velocities. The relatively broader component re- 
flects both the very broad-line cloud ( ^3.11 that happens to 
lie in the same region and the contamination from the much 
brighter emission clouds nearby (EELR-a, b, c). We measure 
the velocities and FWHMs of the two components from the 
[O III] AA4959,5007 lines where we have the best S/N, and 
freeze these parameters when fitting the much fainter [O III] 
A4363 line and the [O ll] doublet. We then use the measured 
line fluxes to derive the electron temperatures and densities 
for individual clouds. 

The electron temperatures of different emission regions ap- 
pear fairly uniform, varying from 12000 to 15000 K. We 
also obtained good measurements of electron densities from 
the relatively isolated clouds (EELR-6, d & e), and we found 
that the density varies from a few tens to a few hundreds 
cm^"^. Lower density gas has higher ionization states (as in- 
dicated by the [O lll]/[0 ll] intensity ratios). The results are 
tabulated in Table 1 . Clearly the pressure is quite different in 
different EELR clouds. 

3.3. Constraints on the ionization mechanism 

Recent Chandra ACIS observations show that 3C 249. 1 has 
extended, highly structured X-ray emission regions, one of 
which follows the structure of the inner optical emission quite 
closely {xa in Fig. nLe/f; region a in Fig. 4 of Stockton et al. 
2006). The X-ray emission of xa cannot be connected with 
X-ray jets frequently seen associated with radio-loud quasars, 
nor can it be due to electron scattering of nuclear emission. It 
could either be thermal bremsstrahlung emission from high- 
speed shocks or X-ray recombination lines from the IC* K 
photoionized gas. If it were the former case, one would expect 
shock-ionized gas near the X-ray emitting region. 

Because of the different shapes of the ionizing-photon spec- 
trum and different temperature regimes, ionization mecha- 
nisms {e.g., shock, "shock + precursor" and pure photoioniza- 
tion) can be distinguished using diagnostic diagrams involv- 
ing line ratios of some optical lines. We measured line fluxes 
from the extracted spectra of various EELR clouds. Figure |6l 
compares the observations with different model grids. Un- 
fortunately, 3C 249.1 EELR clouds (including the region as- 
sociated with X-rays — EELR-d) all fall in the regions where 
the photoionization model and the "shock + precursor" model 
overlap, making it hard to distinguish the two, though the pure 
shock model can be ruled out. We will discuss the X-ray emis- 
sion regions in more detail in §4. 1 . 

We notice that the 3C 249.1 EELR shows systematically 
more Hell emission and lower [Olll] A4363/(A4959h-A5007) 
ratios than 4C 37.43 El, indicating a higher metallicity and 
a lower temperature. But aside from those lines, the spectra 
look remarkably similar to that of 4C 37.43 El . The measured 
line fluxes are tabulated in Table 2, along w ith predictions 
from the two-phase photoionization models of lStockton et alJ 




3.4. Mass of the ionized gas 

The mass of the ionized clouds can be derived from their 
H/3 luminosity and electron density. For the EELR of 
3C 249.1 within the IFU field-of-view, we estimated a total 
H/? flux of 2.9 X 10^^^ erg cm^^ s^^ by assuming a constant 
[Olll] A5007/H/5 ratio of 11 (Fig.|6}. The mass of the H II 
region is then 

47rmp/H/3di 
J^H = r — 1 

where nip is the proton mass, is the luminosity distance, 
aH/3 is the effective recombinat ion coefficient of H/3, and hv 
is the energy of a H/3 photon (Osterbrock 1989). Assuming 
Tie = 100 cm^-^ (§3.2) and case-B recombination at 10"' K, we 
obtain A/h = 6 x lO'^ Mq. A similar calculation for4C 37.43 
El yields a mass of 10^ M©. 

These values are probably underestimates. The mass de- 
rived above assumes a single phase cloud with an average den- 
sity ^100 cm^'^. However, through detailed photoionization 
modeling, Stockton et al. (2002) found that at least two dif- 
ferent density regimes are required to reproduce the observed 
spectrum of 4C 37.43 El : a density -bounded component with 
rie ^ 2 cm^'^and a filling factor / ^ 1 and an ionization- 
bounded component with ^ 500 cm^^^and / ^ 10^^, 
each contributing about one half of the H/3 flux. The spec- 
tra of the EELR clouds of 3C 249.1 look remarkably simi- 
lar to that of 4C 37.43 El (see Table |2jl, so it is appropriate 
to use the same modeling results. Since the [Oil] emission 
comes almost entirely from the high density component, the 
luminosity-weighted electron density ^100 cm^"^ measured 
from [On] AA3726, 3729 doublet (§3.2) favors the densest 
material. This density could be considerably larger than the 
mass-weighted density (i.e., weighting oc rie/) that should be 
used in deriving the mass, therefore the mass derived above 
(6 X 10'' M0) could have been seriously underestimated. For 
example, assuming riei = 2 cm^^ and 71^2 = 150 cm^-^ and 
both components contributing equally to the H/3 luminosity, 
we obtain 1 .5 x 10^ M0 of gas in the low density component 
and only ^2.0x10'' Mq of gas in the high density compo- 
nent. Note that this mass is only an order of magnitude lower 
than the virial mass derived from assuming that the clouds 
are self gravitating and that the measured line widths are en- 
tirely due to motions governed by this gravitational potential 
(i.e., ignoring turbulence or velocity gradients along the line 
of sight): M„„ai - 5i?aVG ~ 3.8 x 10^" R^kpc FWHMIoq 
Mq^. Considering the uncertainties in the simple two-phase 
model (e.g., the possibility of a large neutral fraction), the ac- 
tual masses of individual clouds could potentially be signifi- 
cantly larger than the masses derived from the line fluxes, and 
perhaps the clouds might even be gravitationally bound. This 
could provide an alternative explanation for the existence of 
high-density gas without invoking shocks. 

3.5. Kinetic energy and momentum 

With the knowledge of the mass and kinematics of the 
emission-line clouds, we can estimate their kinetic energy and 
momentum. The bulk kinetic energy of a luminous EELR 
cloud is Eke = Mv"^ /2 = 2.5 x lO^'^MgwIpQ ergs, where 
Afg is the mass of the ionized gas in units of 10^ and 

^ 1 kpc is the average radius of EELR-c measured from the HST image. 
The measurement was corrected for PSF broadening. 
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Fig. 3. — The extraction apertures of different emission-line regions. Baclcground images are the GMOS/IFU [OlII] A5007 radial velocity broad-band channel 
maps. The FOV is slightly different from the one shown in Figs.lTHll because this one is created from the final merged datacube, while the latter is from the 
single exposure which had the best spatial resolution. 



7 
E 

O 

tu 



o 

^ 



0.4 
3 



[0 II] 3726,3729 
/ 



[0 III] 4959,5007 



2 - [Ne V] 342 

; 



1 h 




[Ne III] 3868 



[Ne III] 3967 H? 



I 



I [0 III] 4363 / 




^ 4500 5000 5500 6000 6500 

A (A) 

Fig. 4.— Extracted spectrum of 3C249.1 EELR-fe (refer to Fig.0, with 
important emission lines labeled. 



V500 is the velocity in units of 500 km s^^. The kinetic en- 
ergy of the unresolved kinematic substructures ("turbulent" 
kinetic energy) can be derived from the measured line widths, 
Eturbuient = 6.4 X IQ^^ FWHM§oo crgs. The emission- 
line gas also has a thermal energy ofETH = 2.5 X IQ^^MgTi 
ergs (T4 = T(K)/IQ'^K), which is less than 0.1% of the ki- 
netic energy. Similarly, the total momentum is pkin = Mv 
= 10^*'AfgU5oo dyne s. Whatever the driving source of this 
outflow is, it must have deposited this enormous amount of 
energ y in to the EELR in a short period of time (~10 Myr; 
see ii l4.2> . We w ill try to constrain the power source of this 
outflow in §E21 

4. DISCUSSION 

4.1. The nature of the X-ray emission 

In Fig. we compare the residual Chandra X-ray image 
of 3C 249.1 and the 0'.'35 resolution radio image. This figure 
shows convincingly that one of the X-ray emission region, 
xh, has an intimate relationship with the radio jet. The X-ray 
emission is coincident with a region protruding perpendicular 
to the jet direction, a feature that has been n oted by several au- 
thors {e.g., designated XI in Gilbert et al. 2004). Since there 
is a large offset (~1'.'4 Ki 6.4 kpc) between xh and the nearest 
hotspot (N2, foUowing Rjilbert et al1l2004 . it might not sim- 
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cussed in §3.1 and the contamination from the much brighter emission clouds 
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dashed curves and dotted curves show the "broad" and "narrow" component, 
respectively. The fit to the [Oil] doublet have been constrained so that the 
relative velocity between the "broad" and "narrow" components are frozen to 
that determined from the fit to the [O III] lines. 



ply be synchrotron radiation from the jet, but possibly a result 
of a gas cloud being shocked by the jet fluid. 

The other X-ray emission region, xa, is covered by the IFU 
FOV. Gas kinematics, line-ratio diagnostics and its morphol- 
ogy all fail to provide convincing evidence of shocks, so we 
can not rule out the possibility that the X-ray emission is actu- 
ally a result of the same recombination process as what causes 
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Fig. 6. — Optical line-ratio diagrams are oveiplotted with model grids. The dusty radiation pressure-dominated photoionization model (Groves et al. 2004), 
shock-only and "shock + precursor" model (Dopita & Sutherland 1996) grids are plotted in grey, blue and green, respectively. The photoionization grids assume 
a power-law photoionizing continuum with index a = —1.4 and show a range of ionization parameters (Logio(U)) for metallicities of 0.5 Zq, 1.0 Zq and 
2.0 Zq . The shock + precursor models assume equal contribution to the H/3 line from the shocked gas and the precursor A range of shock velocities (Vg [km 
s~^]) and magnetic pai'ameters (B/n^^^ [fi G cm~^/^] ) are covered. Line ratio measurements fro m four major EELR c louds of 3C 249.1 are shown in red, 
and the blue point shows the measurement from the Keck II/LRIS long-slit spectrum of 4C 37.43 El iStockton et alJ2003) . Arrows are reddening vectors of 3 
Ay, assuming the standard Galactic reddening law and Rv = 3.1. 
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Fig. 7.— PSF-subtracted Chandra ACIS image of 3C 249.1 overlaid with 
contours of a VLA image at 5 GHz and 0'.'35 resolution ( Gilbert et al. 2004|). 
Contours at 0.5 X (1, 2, 8) mjy beam^^. The two X-ray emission regions (xa 
& xb) and a radio hotspot (N2) are labeled. The peak of the radio continuum 
has been registered to that of the X-ray emission, which is indicated by a 
cross, before a scaled synthetic PSF was subtracted from the X-ray image 
(refer to iStockton et all2006l for more details on the PSF subtraction). 



the optical emission. If the X-ray is indeed from recom- 
bination hnes, then the other high surface-brightness clouds 
(EELR-a, b, c) should also be bright in X-ray; but it is im- 
possible to detect the emission from the Chandra data, since 
at such small radial distances their X-ray emission, if any, is 
completely overwhelmed by the X-rays from the quasar. 

4.2. Constraints on the power source of the outflow 

The dynamical timescale of the EELR can be estimated 
from their radial velocities and distances to the nuclei, r^yn ~ 
10^Z)5fepc^^500 yr- Assuming an average velocity of 500 km 
s~^ is appropriate here, because the initial velocity of a galac- 
tic wind is expected to be on the level of 1000 km s^^ and it 
has since been slowed to the observed a few hundreds km s^^ 
due to mass entraining and working against gravity. At the 
redshift of 3C 249.1, the angular scale is 4.6 kpc arcsec^^. 



so 5 kpc is about the projected distances of the two bright- 
est EELR clouds to the nucleus. The time scale is consistent 
with the age of typical luminous radio lobes, ^ 10^ — 10^ 
yr ( Blundell & Rawlinas 2000). The rate of mass outflow 
can be estimated assuming it is constant over the dynamical 
age of the ionized gas, M ~ M/rdyn = 100^^9^5fcpc^500 
M0 yr^^. Similarly the injection rate of kinetic energy is 



Eke 

-,-1 



'5fcpc^^500 ergs 

s^"^, and the input momentum rate is pkin — Pkin/Tdyn — 
3 X IO^^MqV5oqD^^p^V5oo dyne. We emphasize that the rates 
can only be considered as lower limits to the actual wind en- 
ergetics since the majority mass of the wind flow may not be 
detectable in the optical ( Veilleux et al. 2005) and a substan- 
tial amount of wind-entrained ISM may be shielded from the 
ionizing flux of the quasar. 

If the outflow is driven by a starburst, a star formation rate 
(SFR) of ~10 M0 yr^^ is required to deposit enough kinetic 
energy into the clouds, or ^60 Mq yr^^ to inject e nough mo- 
mentu m into the ionized gas (Equations 2-3 in V eilleux et alJ 
120051 Since a significant fraction of the originally injected 
kinetic energy may have been lost to radiation and counter- 
acting the gravitational potential, the kinetic energy is not 
conserved. But the wind may still conserve momentum, so 
the SFR of ~60 Mq yr^^ is probably more realistic. The 
SFR scales to an infra-red (IR) luminosity of L/j7(8-1000/^m) 
« S.SxIO^SFRLq = 3.5x10"Lq (j^nnicutt 1998). Is this 
luminosity consistent with mid-IR observations? 3C 249.1 
was d etected in 12, 25 and 60 ^m by IRA S (|Sanders et aQ 
[T989I) . in 70 fim by MIPS (Shi et al. 2005), and in 100 jim 
by ISO ((Haas et al, 2003). By using a model of two black- 
body components to fit the mid-IR spectral energy distribu- 
tion (SED), we obtain a total IR luminosity of L/7j(8-1000 
^m) ~ 8.6 X 10" Lq. Thus more than 40% of the total IR 
luminosity must be powered by star formation to explain the 
input momentum rate. However, the SED peaks near 12 /im, 
so the IR emission is dominated by warm dust (T 100 — 200 
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K) presumably directly heated by the quasar JSanders et alJ 
[1289). From the SED fitting, less than 23% of the total lumi- 
nosity could have come from a "starburst" component with a 
characteristic dust temperature T < 65 K. Apparently the IR 
data argues against an on-going starburst as strong as ~60 M0 
yr^^, though it doesn't formally rule out the possibility that 
there was an energetic starburst which had caused the outflow 
but ceased not too long ago (< Tdyn ~ 10 Myr). 

Note also that the mass outflow rate is almost an order of 
magnitude larger than the predicted mass injection rate from 
supernovae (SNe) for a SFR-60 M© yr'^ (Af,„j ~ 16 M© 
yr^^), which implies that the wind has entrained considerable 
amount of gas from the ISM. This is in contrary to the low 
efficiency for transporting mass out of the galaxy predicted 
by numerical simulations o f starburst-driven galactic winds 
JStrickland & Stevensl2000l) . 

Now we consider the possibiUty that the outflow is di- 
rectly driven by the quasar of 3C 249.1. Radiation from the 
quasar can couple to the surrounding gas via various pro- 
cesses such as electric scattering, photoionization, etc. The 
input momentum rate from radiation pressure is p = 1.3 x 
10^^(77/0. l)Macc dyne s, where 77 is the radiative efficiency 
and Mace is the mass accretion rate of the black hole in units 
of M0 yr~^. Hence a mass accretion rate of 2.5 Mq yr^^ is 
sufficient to inject enough momentum to the clouds in 10^ yr 



This accretion rate is in agreement with the one implied by its 
bolometric luminosity. Mace = 8(77/0.1)^^ Mq yr^^ (Lboi = 
8 X 10^2 Lq; Sanders et al. 1989). 

In summary, although we cannot totally exclude a starburst 
origin, the available evidence favors the quasar as the power 
source of the outflow. If this is indeed the case, the presence 
of a lu minous EELR may tell us so mething about the quasar 
itself. Stockton & Mac Kentv* ('1987) found that roughly one- 
third to one-half of all powerful, steep-spectrum radio-loud 
quasars show luminous EELRs, but there are clearly some 
for which any extended emission is quite weak. These lat- 
ter cases may be those for which the quasar ignition occurred 
sufficiently long ago (> a few lO*^ years) that the expelled 
gas has dissipated to the point that its surface brightness falls 
below our detection limit. 
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TABLE 1 

Electron Density and Temperature of the EELR Clouds. 





V 


FWHM 








Region 


(kms-i) 


(kms-i) 


(K) 


(cm-3) 


[Olll]/[Oll] 


3C 249 B 


-90 


300 


12000+ 800 


160 + 100 


5.6 


3C 249 D 


+140 


300 


12100 ± 1300 


300 ± 140 


4.2 


3C249E 


+300 


120 


15000 ± 1500 


30 ± 100 


7.2 



TABLE 2 

Observed Line-Flux Ratios for 3C 249.1 EELR Clouds and 4C 37.43 El and Modeled Line-Flux Ratios 



Line Flux Ratios" 


[Ne V] A3426 


[0 II] AA3726, 9^ 


[Ne III] A3869 


[Olll]A4363 


He II A4686 


[Olll] A5007 


3C249 A 


0.44 + 0.15 


2.39 + 0.17 


1.00 + 0.12 


0.16 + 0.07 


0.29 + 0.08 


10.80 + 0.50 


3C249B 


0.60 + 0.15 


2.22 + 0.14 


0.98 + 0.11 


0.17 + 0.07 


0.32 + 0.07 


12.47 + 0.46 


3C249C 


0.30 + 0.13 


2.52 + 0.19 


0.89 + 0.12 


0.14 + 0.06 


0.27 + 0.06 


10.89 + 0.46 


3C249D 


0.63 + 0.40 


2.47 + 0.24 


0.71 + 0.13 


0.14 + 0.09 


0.28 + 0.08 


10.38 + 0.54 


4C37E1 


0.34 + 0.04 


2.86 + 0.07 


0.95 + 0.05 


0.18 + 0.05 


0.20 + 0.03 


10.50 + 0.22 


Modell' 


0.34 


1.77 


0.85 


0.19 


0.20 


10.31 


Model2 


0.30 


2.33 


0.81 


0.17 


0.20 


9.35 



"All line fluxes are given as ratios to the H/3 flux 
''Total flux in the [O II] doublet 

■^Models 1 and 2 have 25% and 33% of the H^S flux coming from the 500 cm"^ high-density component, see Stockton et al. 2002 (§3.4.2) for details 



